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A NOVEL METHOD OF GAMMA RAY EVENT SELECTION 
Field Of The Invention 

[0001] The present invention relates to an improved 

methodology of discriminating between desired and 
undesired readings or events measured by radiation 
detectors. Specifically, the method is intended for use 
in radiation detectors having non-Gaussian energy 
response functions. 

Background Of The Invention 

[0002] Gamma ray cameras are relatively well known 
devices used to detect gamma ray emissions from 
radioactive decay. A known gamma ray camera described in 
United States Patent No. 3,011,057 for RADIATION IMAGE 
DEVICE, hereby incorporated by reference, uses a single 
sodium iodide ("Nal") scintillation detector crystal to 
detect gamma ray emissions. The detector crystal is 
positioned to receive a portion of the gamma ray 
emissions from the radioactive decay. When a gamma ray 
strikes and is absorbed in the detector crystal, the 
energy of the gamma ray is converted into a large number 
of scintillation photons that emanate from the point of 
the gamma ray's absorption in the detector. A 
photomultiplier tube, optically coupled to the detector 
crystal, detects a fraction of these scintillation 
photons and produces an electronic signal that is 
proportional to the number of detected incident 
scintillation photons. The gamma ray camera typically 
has several photomultiplier tubes placed in different 
positions, with the signals from the different 
photomultiplier tubes being combined to provide an 
indication of the positions and energies of detected 
gamma rays. 



[0003] Gamma ray cameras are frequently used in 
nuclear medical imaging to record pictures of a radiation 
field emanating from a subject's body. The gamma rays 
originate from a decay of a radioactive tracer that has 
been introduced into the subject's body. The radioactive 
tracer, such as 99m Tc, is a pharmaceutical compound to 
which a gamma ray emitting nuclide has been attached and 
which undergoes some physiological process of interest 
after introduction into the body. For example, the 
tracer may accumulate in areas of high blood flow, 
thereby pinpointing areas of physiological activity. 
[0004] Only substantially collimated gamma rays reach 
the gamma ray camera because a collimator usually is 
placed between the source of radiation and the 
scintillator. Of these collimated gamma rays, the gamma 
ray camera only detects the small fraction of gamma rays 
that impinge a detector, such as the above-described Nal 
detector crystal and photomult iplier tube assembly. From 
these detected gamma rays, the gamma ray camera typically 
selects a desired sample of gamma rays based upon their 
measured energy according to a process described in 
greater detail below. The gamma ray camera then records 
in an image memory the number of gamma events detected at 
each of a number of spatial locations corresponding to 
regions in the patient's body from which the gamma rays 
emanated. The data in the image memory is then used to 
form an image that corresponds to the distribution of the 
detected gamma events, thereby providing an image of the 
organ, tissue or body region of interest. 

[0005] The quality of the gamma ray camera's image of 

the distribution of the radioactive tracer is dependent 
on the sensitivity of the detectors in the camera. As 
the sensitivity of the detector to gamma rays is 
increased, more gamma rays are detected and incorporated 
into the image, because the overall detection system 
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sensitivity, combined with the activity of radioactive 
tracer presence in the body, determines the number of 
detected gamma ray events. Creating an image from a 
small set of gamma ray events may result in an inaccurate 
image because the small sample may not accurately 
represent the true distribution of the radioactive tracer 
within the subject's body. In particular, a reduction in 
the number of detected events, or "counts," results in 
increased fractional statistical uncertainty because the 
emission and subsequent detection of gamma-rays are both 
random, stochastic processes. As a result, for any 
single gamma ray detector, the standard deviation in the 
number of counts in that detector increases, as a 
fraction or as a percentage, as the number of counts 
decreases. Therefore, the image would be more noisy in 
terms of image quality as a result of decreasing the 
total number of counts contributing to the image. 
[0006] However, increasing the total number of counts 
by the inclusion of false events may degrade the image 
produced by the camera. One source of false events is 
the scattering of some of the gamma rays between the 
point of emission and exiting the patient's body. These 
scattered gamma rays do not provide an accurate 
indication of the distribution of the radioactive tracer 
within the patient ' s body because the scattering changes 
the gamma ray's direction and the measured energy. As a 
result, the inclusion of too many scattered gamma rays in 
the image memory causes the gamma ray camera to produce 
an inaccurate image with poor contrast. 

[0007] Gamma ray detectors, therefore, seek to strike 

a balance between scattered and unscattered gamma rays 
without a loss of sensitivity. Typically, this result is 
accomplished by evaluating the energy levels of the 
detected gamma ray events and selecting those gamma ray 
events that fall within a preset "window" or range of 
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acceptable energy levels. This technique is based on the 
assumption that gamma rays of abnormally low or high 
energy levels have been altered prior to detection or did 
not originate from the radioactive tracer of interest 
and, accordingly, do not provide reliable indications of 
the rays' direction of origination. The unscattered 
gamma ray events tend to fall within the window, whereas 
the scattered gamma ray events tend to fall outside the 
window. While a small number of scattered gamma rays 
would also fall within this energy window, these detected 
scattered gamma rays would be greatly outnumbered by the 
unscattered events. The gamma ray camera may vary the 
energy window as needed to achieve a desired result; for 
example, a narrower energy window may continue to accept 
a large number of unscattered gamma ray events while 
excluding more scattered events. 

[0008] Frequently, the energy resolution of the gamma 

ray detector, such as the above-described Nal detector, 
can be characterized by a Gaussian distribution. The 
energy resolution of the gamma ray detector is then 
described by the full-width at half-maximum ( " FWHM" ) of 
the Gaussian shaped photopeak. 

[0009] As explained above, gamma ray cameras generally 
use energy windows to discriminate between scattered 
gamma rays and gamma rays that reached the detector 
without interaction. The use of the energy window that 
defines a range of acceptable energy events is fairly 
effective in gamma ray cameras having Gaussian energy 
response functions because a large percentage of the 
desired gamma ray events are concentrated within a small 
range of energy levels. For instance, with a 
monoenergetic gamma ray emitting isotope in the absence 
of any scattered gamma rays, an energy window of 
approximately twice the width of the FWHM centered over 
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the photopeak would select on the order of 99% of the 
total unscattered events. 

[0010] However, the use of an energy window has the 
unavoidable disadvantage of excluding some of the desired 
gamma ray events while including some of the undesired 
gamma ray events, regardless of the width of the energy 
window. It is therefore the goal of the present 
invention to present an improved methodology to 
differentiate between the desired and undesired gamma ray 
events, without decreasing the sensitivity of the 
detector. 

[0011] U.S. Patent No. 5,561,297 for SCATTER 

CORRECTING GAMMA CAMERA, incorporated by reference herein 
in its entirety, discloses a method that improves the 
quality of an image by reducing the inclusion of 
scattered gamma rays into an image. The method measures 
the response of the imaging detector to substantially 
unscattered gamma rays, fitting a calibration function to 
the detector response to the unscattered gamma rays for 
each spatial location. The method then fits the response 
of the imaging detector to a field of both scattered and 
unscattered gamma rays with a combination of the measured 
responses to the unscattered gamma rays and an estimated 
shape for the spectrum of scattered gamma rays. The 
response function is a Gaussian function and the centroid 
energy and standard deviation of the Gaussian function 
are the spatially varying parameters that are determined 
by fitting the unscattered gamma ray spectrum for each 
spatial location. 

[0012] The disadvantages of using an energy window are 
especially evident with gamma ray detectors having non- 
Gaussian energy response functions. In detectors with an 
energy response function that is not Gaussian shaped, a 
significant number of the desired gamma ray events occur 
away from the photopeak of the energy response function. 
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The use of a narrow energy window would exclude these 
desired gamma ray events, whereas the use of a wide 
energy window to include these off-peak values would 
include too many of the undesired, scattered gamma ray 
events . 

[0013] One type of gamma ray detector that has a non- 
Gaussian shaped energy response function is a solid 
state, pixelated detector made from the room temperature 
semiconductor material cadmium zinc telluride ("CZT"). 
For example, the measured energy response function for 
one CZT detector has a photopeak in the range of 5% FWHM 
and a significant tailing toward the lower energy levels. 
Therefore, a large number of the desired gamma ray events 
detected with this CZT detector occur away from the 
photopeak and would be wrongly excluded with the use of a 
narrow energy window. 

[0014] FIG. 1 illustrates a typical energy spectrum of 

a single pixel of a pixelated CZT detector exposed to 
substantially unscattered 140 keV gamma rays. In FIG. 1, 
the peak centered about 140 keV represents the gamma rays 
that have been absorbed substantially within the center 
portion of the single pixel, and the distribution of 
signal amplitudes of these events is approximately 
Gaussian. However, a significant number of gamma rays 
are also detected in the portion of the energy response 
spectrum that tails toward the lower energies. This 
phenomenon is caused, in part, by gamma ray absorption 
events that do not confine all charge creation to within 
a single pixel and by non-ideal charge collection. 
Because the illustrated response function represents the 
distribution of measured signals from only a single 
pixel, charge that is lost from the pixel and shared with 
adjacent pixels is not included in the response function. 
As a result, gamma ray absorption events in which the 
charge collection is incomplete due to charge sharing 
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with other pixels are lost from the peak region and 
contribute to the low energy tailing. 
[0015] The use of a simple energy window to 
discriminate unscattered from scattered gamma-rays for 
detectors with non-Gaussian response functions, such as 
the CZT detector, is problematic. The use of a simple 
energy window for the detector pixel shown in FIG. 1 
would reduce the overall sensitivity of the detector if 
the window width were made narrow (about two times the 
FWHM width of the peak) . A narrow window would encompass 
the peak but not the unscattered events that contribute 
to the tailing toward lower energies, so a large number 
of the desired gamma ray events would be excluded, 
reducing the sensitivity of the detector. Conversely, 
using a very wide energy window to include the peak 
region and the tailing toward lower energies would also 
include a high fraction of scattered gamma rays in the 
image, and thus would degrade the image quality. 
[0016] It is therefore a further goal of the present 

invention to provide a reliable methodology for 
estimating unscattered gamma rays detected by a gamma ray 
detector having a non-Gaussian energy response function. 

Summary Of The Invention 
[0017] These and other needs are addressed in the 
present invention through a novel method for 
discriminating between scattered and unscattered detected 
radiation events, so that the sensitivity of the 
radiation detector may be preserved without causing the 
inclusion of a large number of undesired radiation 
events. This method includes the steps of: (1) 
empirically determining the energy response function of 
the radiation detector; and (2) then using the 
empirically determined energy response function during 
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the actual detection process to differentiate between 
desired and undesired radiation events. 

[0018] In a preferred embodiment, the energy response 

function for each of the detectors is determined by 
simulating the condition for the subsequent, actual 
measurement. During this process, the detector is 
illuminated with a uniform field of radiation to be 
measured. The detector should be exposed to a large 
number of radiation photons, and the number of undesired 
rays is minimized. For example, the detector will be 
exposed to only direct, unscattered radiation. This 
simulation produces an energy response function for the 
detector under conditions that approximate the actual 
measurement . 

[0019] During the actual measurement, the detector 

produces a second energy response spectrum that 
represents both the desired and undesired events. Then, 
a least squares estimate of the number of desired, 
unscattered events is produced by taking the dot product 
of (1) the spectrum acquired during actual measurement 
and (2) a weighting vector determined during the 
calibration step. 

Brief Description Of The Drawings 
[0020] These and other features and advantages of the 
present invention are more fully described in the 
following attached drawings in which corresponding 
elements are designated by like reference numbers: 
[0021] FIG . 1 (Prior Art) illustrates an energy 

spectrum of one type of detector having a non-Gaussian 
energy response function; 

[0022] FIGS. 2-12 illustrate energy spectrums from a 

hypothetical example illustrating the method of FIGS. 
13A-C; and 
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[0023] FIGS. 13A-C are flow charts of the method for 

discriminating between detected radiation events in 
accordance with an embodiment of the present invention. 

Detailed Description Of The Preferred Embodiment 
[0024] As illustrated in FIG. 13A, the present 
invention is a method of replacing the simple energy 
window as a means for discriminating between scattered 
and unscattered gamma rays in a gamma ray camera. The 
method has two steps: 1) a calibration step 100 to obtain 
the energy response function of the detector; and 2) an 
acquisition and correction step 200 when radiation is 
actually detected and then selected for use according to 
information determined during the calibration step 100. 
Details of the calibration step 100 are shown in FIG. 
13B, and details of the acguisition step are shown in 
FIG. 13C. 

[0025] The response function of each individual pixel 
in an imaging detector is measured, giving a result 
similar to the energy response function in FIG. 1, at 
step 150 of FIG. 13B. Measurement of the response 
function is accomplished by illuminating the detector 
with a uniform stream of radiation of the energy or 
energies of interest, step 110. The uniform flood 
ideally comprises only unscattered gamma rays, so the 
number of scattered radiation events should be minimized, 
step 120. During the calibration step, 110, the total 
number of counts acquired by each detector should be 
large, so that the shape of the response function is 
accurately measured and not dominated by random counting 
statistics, step 130. 

[0026] It is also important that the conditions of the 

illumination during the calibration be substantially 
identical to conditions encountered during actual 
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measurement, step 140, such as the conditions encountered 
during the clinical imaging of a patient's body when a 
large fraction of scattered gamma rays will be present. 
For example, the response function is generally 
significantly different for a pixelated detector 
illuminated uniformly over its entire surface as compared 
to illumination confined to the central region of each 
pixel detector. Therefore, during the measurement of the 
response function in the calibration step 100, the 
radiation collimation must match the collimation to be 
received by the detector during the actual measurement in 
step 200, so that the response functions will remain 
constant . 

[0027] Alternatively, the detector may use previously 
measured radiation samples and use these results in 
calibration step 100 to estimate the energy response 
function for the detector. 

[0028] Where the detector is pixelated, such as the 

above described cadmium zinc telluride ("CZT") detector, 
a separate energy response function is determined for 
each pixel during the calibration step 100. 
[0029] During the acquisition step 200 of the method, 
as shown in FIG. 13C, the detector makes actual 
measurements, step 210, and forms energy spectra for 
these measurements, step 220. The radiation measurements 
include both desired and undesired events. During 
acquisition, detected events are stored in a memory of 
the gamma ray camera. The acquisition can either acquire 
a large list mode file, where all events are stored 
sequentially with the address of a detector or event 
location and the energy or amplitude of that event, or 
individual energy spectra can be accumulated for each 
pixel detector separately. 
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[0030] The spectrum corresponding to the response of 
the detector to unscattered gammas of energy, E, is 
defined as the vector, R E . The vector length is the 
number of energy bins in the spectrum. Similarly, the 
spectrum corresponding to the response of the detector 
pixel to scattered gammas from a source of energy, E, is 
defined as the vector, S E . Using this convention, the 
number of events or counts detected in each energy bin 
within a single detector is represented as a vector C E . 
The spectrum recorded during an acquisition is then given 
by Equation 1: 



In Equation 1, c us is the number of unscattered counts and 
c s is the number of scattered counts acquired, and 
accordingly, c us and c s are scalar quantities. 
[0031] Equation 1 may be rewritten as Equation 2, 




where F is given by Equation 3, 
R ? S ? 

F= 2 2 (3) 
R S 

\ n n) 

Solving for scalar coefficients yields Equation 4, 

= ( F T F) _1 F T C E (4) 
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In this solution, the matrix (F T F) """F 1 consists of two 
rows. The first row is a column vector V 1 that represents 
the weighting factors that correspond to the unscattered 
component in a spectrum. Therefore, the dot product of 
the first row vector V x with a column vector spectrum C E 
equals the least squares estimate of the magnitude of 
unscattered component in the spectrum, which is c us . The 
dot product of the acquired spectrum C E in step 120 with a 
weighting vector V x determined during the calibration 
gives a least squares estimate of the number of 
unscattered gammas in the acquired spectrum, step 230. 
This number c us is exactly what is desired for determining 
how many of the measured events to include into the 
image . 

Example 1 

[0032] For the purpose of development and 

demonstration of the method of this invention, a Monte 
Carlo simulation was used to generate spectra of 
scattered gamma rays emanating from a sphere of water. 
[0033] FIGS. 2-6 show the components and composite 

spectrum emanating from a 10 cm diameter water sphere 
with a 140 keV source of gammas at the center. Estimates 
of the measured spectra for unscattered and scattered for 
the CZT detector pixel were derived by convolving the 
Monte Carlo results with the detector response function 
measured at 140 keV. Specifically, FIG. 2 illustrates 
the energy response spectrum of single Compton scattered 
gamma rays; FIG. 3 illustrates the energy response 
spectrum of twice Compton scattered gamma rays; FIG. 4 
illustrates the sum spectrum from FIGS 2-3 for the 
scattered gamma rays, normalized to 0.78; FIG. 5 
illustrates a monoenergetic source spectrum normalized to 
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0.22; and FIG. 6 illustrates a simulated total spectrum 
from FIGS 2-5. 

[0034] FIGS. 7-9 respectively show the unscattered, 

scattered, and composite spectra expected for the CZT 
detector, assuming an arbitrary mixture of 70% scatter 
and 30% unscattered gamma rays incident on the detector. 
Testing methods of determining the number of unscattered 
gamma detections were done by combining arbitrary 
fractions of scattered and unscattered gammas, and then 
randomly generating a detected spectrum for a fixed 
number of total counts. Specifically, FIG. 7 illustrates 
an unscattered spectrum in a CZT detector normalized to 
30%; FIG. 8 illustrates a scattered spectrum in a 
detector normalized to 70%; and FIG. 9 illustrates the 
composite spectrum for the CZT detector from FIGS 7-8. 
Figure 10 then shows an example spectrum containing 100 
counts, randomly selected from the composite spectrum of 
FIG. 9. 

[0035] Three test spectra were then constructed by 
adding together different proportions of scattered and 
unscattered response. The three choices were: 1) 50% 
scatter/50% unscattered, 2) 70% scatter/30% unscattered, 
and 3) 90% scatter/10% unscattered. For each of these 
cases, four trials were done where the total numbers of 
counts in a trial were fixed at 500, 100, 50, and 10 
counts. For each of these trials, 100 random 
distributions were generated using a Poisson random 
number generator. Each of these random spectra was 
similar to Fig. 10. 

[0036] A comparison was made between the method of 
this application and the conventional method of using a 
simple energy window to determine the number of 
unscattered events to include in the image. For each 
trial, each of the 100 random spectra was used to compute 
the counts that a pixel would record in the image. A 
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fixed window of 10% width (14 keV total, centered at 140 
keV) was simulated by simply summing all of the random 
counts that fell within the window. A 10% window was 
chosen because the peak in the response function is 
approximately Gaussian with a 5% FWHM. The method of 
energy discrimination of this application was simulated 
by simply computing the dot product of the unscattered 
component vector with the random counts spectra. The mean 
number of counts and the standard deviation of each trial 
were computed, for both the window and the method of this 
invention using the 100 random samples for each method. 
The results are summarized in Table 1. 



TABLE 1. COMPARISON OF FIT METHOD TO SIMPLE ENERGY WINDOW 



[0037] The results summarized in Table 1 clearly show 
that the method of this invention is far superior to an 
energy window for predicting the mean number of 
unscattered detection events for a detector of the 
response function of the subject CZT detector. Although 
the estimates using the window method have a smaller 
variance about their predicted mean, their predicted mean 
is significantly in error in all cases. 

[0038] FIG. 11 illustrates the vector corresponding to 

the weights of the unscattered component from the 
measurement; and FIG. 12 illustrates the vector 
corresponding to the weights of the scattered component 
from the measurement. The fit or weighting vectors 
illustrated in Figs. 11 and 12 were derived using the 70% 
scatter/30% unscattered spectrum. 

[0039] It is important to note that the predictive 

accuracy of the unscattered counts is almost as good for 
the 50% scattered/50% unscattered and 90% scattered/10% 
unscattered cases, as for the 70% scattered/30% 
unscattered. This fact suggests that the method is 
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tolerant over a large range of variability in the scatter 
to unscattered components. 

[0040] The invention having been described, it will be 
apparent to those skilled in the art that the same may be 
varied in many ways without departing from the spirit and 
scope of the invention. Any and all such modifications 
are intended to be included within the scope of the 
following claims. 
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